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The invariant mass spectra of e + e~ pairs produced in 12-GeV proton-induced nuclear reactions 
are measured at the KEK Proton Synchrotron. On the low-mass side of the w meson peak, a 
significant enhancement over the known hadronic sources has been observed. The mass spectra, 
including the excess, are well reproduced by a model that takes into account the density dependence 
of the vector meson mass modification, as theoretically predicted. 
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It is well established that most of hadron masses are 
generated due to the spontaneous breaking of the chiral 
symmetry, which is the crucial aspect of the strong in- 
teraction. The modification of hadron masses and decay 
widths in hot and/or dense matter is theoretically pre- 
dicted as a consequence of the restoration of the broken 
symmetry. Experimental observations of such phenom- 
ena have become one of the most interesting topics in 
hadron physics today. 

As for theoretical studies, many works have been per- 
formed and they are summarized in 0,0- Specifically, 
our experiment was motivated by two related studies 
concerning dense matter. Using an effective chiral La- 
grangian, Brown and Rho proposed an in-medium scaling 
law, which predicted a decrease of the vector meson mass 
at 20% at the normal nuclear density po 0] ■ Hatsuda 
and Lee calculated the density dependence of the mass of 
vector mesons based on the QCD sum rule to reach the 
conclusion that the mass is approximately linear to the 
density in < p < 2po, and the mass of p and u> mesons 
decreases at about 16% at po 00]. 

The experiment E325 has been performed at the KEK 
12-GeV Proton Synchrotron to measure the invariant 
mass spectra of p,co,4> — > e + e~ and <j> —>K + K~ decay 
modes simultaneously. Our main goal is to detect the 
modification of the spectral shape of vector mesons in 
nuclear media. In our earlier publication, we showed the 
significant enhancement below the oj peak in the invariant 
mass spectrum of e + e~ pairs for a copper target, which 
is attributed to the mass-modification effect in nuclear 
matter 0]. This was the first dilepton measurement of 



the modification of the vector meson mass in nucleus. 

It should be noted that our earlier and present re- 
sults are strongly related to other few experimental obser- 
vations. The CERES/NA45 collaboration observed the 
low-mass electron-pair enhancement in Pb-Au collisions 
at 158 A GeV @. The STAR collaboration reported 
the shift of p meson peak in the ir + n~ channel at RHIC 
in Au+Au and p+p collisions at y/s = 200 GeV/c 2 0]. 
Their results may connote the chiral symmetry restora- 
tion in hot matter. The TAGX collaboration reported 
an in-medium modification of the p invariant mass distri- 
bution via 3 He, 12 C(7, 7r + 7r _ )X reactions 10]. Recently, 
the CBELSA/TAPS collaboration has reported that the 
mass modification of w meson in the ir°j channel flU . 

Although these diverse experimental signatures, in- 
cluding ours, are not ample, they could be related to 
a critical behavior of hadrons in hot and/or dense mat- 
ter. Compared to other experiments, our experiment is 
unique, since we are able to measure the e + e~ mass 
shapes directly with a high mass resolution and high 
statistics, placing a special emphasis on vector mesons 
that decay inside a nucleus. Needless to say, lepton pairs 
are almost free from final state interactions with nuclear 
matter. 

The results presented in this paper are based on most 
of the data that we have acquired. The statistics have 
improved by about 50-times compared to those given in a 
previous paper 0] ■ We will show a comparison of the data 
to the theoretical prediction using a toy model dealing 
with the density dependence of the vector meson masses. 

A detailed description of the E325 spectrometer can 



2 



be found elsewhere |l2J. The detector elements, which 
are relevant for this analysis, are briefly discussed be- 
low. The spectrometer was designed to detect decays 
of slowly moving vector mesons, which have a larger 
probability to decay inside a target nucleus. In 2002, 
the primary proton beam with a typical intensity of 
7xl0 8 Hz was delivered to one carbon and four cop- 
per disk targets, which were aligned in-line along the 
beam axis. The interaction length of each copper target 
is 0.054%(73 mg/cm 2 , 0.57% in radiation length), and 
that of the carbon target is 0.21%(184 mg/cm 2 , 0.43% 
in radiation length). The tracking system, consists of 
three drift chambers, gives the momentum resolution of 
a p /p = v /(l-37%-p) 2 + 0.41% 2 [13. 

The mass resolution and mass scale were examined for 
the observed peaks of A — > p7r~ and K° — > it + it~ de- 
cays. The observed resolutions and centroids of these 
resonances are consistent with the expectations given by 
a detailed detector simulation using Geant4 hj. The 
effects caused by multiple scattering, energy losses, the 
chamber resolution, and a miss-alignment of the track- 
ing devices were minutely inspected. Please note that 
the mass spectra presented in this manuscript is not cor- 
rected for the mass scale which is 3.7±0.8 MeV/c 2 lower 
than the real value at 0.8GeV/c 2 (i.e. to peak). At 
the low-mass side of the resonances, a long tail arises 
due to the energy loss of electrons, mainly caused by 
bremsstrahlung. Such a tail is estimated to be negligible 
compared to the excess that we have been observing. The 
mass resolutions for lu — > e + e~ and <f> — ► e + e~ decays are 
estimated to be 8.0 and 10.7 MeV/c 2 , respectively. These 
values are consistent with the data when we neglected the 
excess parts of the mass spectra. 

Figure n shows the invariant mass spectra of the e + e~ 
pairs using all of the data taken in 2002. We have re- 
quired each track of an e + e~ pair to come into each of the 
two arms; therefore, the low-mass region of the spectra is 
largely suppressed. The invariant mass spectra was fitted 
with the combinatorial background and known hadronic 
sources: p — ► e + e~, lu — ► e + e~, <j> — ► e + e~, r) — ► e + e~j, 
and lu — > e + e~7r°. The combinatorial background was 
evaluated by the event-mixing method. The relativistic 
Breit-Wigner distribution was used to obtain the spectral 
shapes of resonances. The mass resolution and the de- 
tector effect of our spectrometer were taken into account 
through the detector simulation described before. The 
kinematical distributions of mesons have been obtained 
by the nuclear cascade code JAM Q , which is in a good 
agreement with the real data. The relative abundances 
of these components were determined by the fitting, ex- 
cept for the ratio of lu — ► e + e~7r° to lu — > e + e~ decay 
which was fixed to their branching ratios, 59/6.95 given 
by the PDG 0] . The fit results are plotted with the solid 
lines in FigFJand summarized in Table |U The obtained 
X 2 /d.o.f. arc 161/140 and 154/140 for the C and Cu tar- 
gets, respectively. The region from 0.6 to 0.76 GeV/c 2 
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FIG. 1: Invariant mass spectra of e + e~ for the (a) C and (b) 
Cu targets. The solid lines are the best-fit results, which is 
the sum of the known hadronic decays, lu — > e + e~ (dashed), 
4> — » e + e~ (thick dash-dotted), r\ — > e + e~j (dash-dotted), 
and cu — > e + e~7r° (dotted) together with the combinatorial 
background (long-dashed), p — > e + e~ is not visible (see text). 



TABLE I: Signal yields in the acceptance together with 
X 2 /d.o.f. for the C and Cu targets, obtained by the fit ex- 
cluding the mass range of 0.6 to 0.76GeV/c 2 . The values for 
p are expressed as an upper limit at 95% C.L. 



77 Dalitz 



xVdof 



C 1012±112 398±42 3644±92 (112) 1461±131 161/140 
Cu 1249±126 547±45 3346±91 (169) 1341±136 154/140 

TABLE II: Same as Table|l]but corresponds to the fit includ- 
ing the excess region. 



r\ Dalitz 



C 1025±116 347±40 3552±98 1468±242 893±169 366/162 
Cu 1248±132 505±44 3209±98 1405±257 798±175 295/162 



was excluded from the fit, because the fit including this 
re gion resulted in failure at C.L. 99.9% as listed in Table 

MM- 

A significant excess can be seen on the low-mass side 
of the lu peak, whereas the high-mass tail of the lu can 
be reproduced with the expected shapes. The number 
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FIG. 2: Invariant mass spectra of e + e~. The combinatorial background and the shapes of 77 — ► e + e~j and u) —* e + e~n 
were subtracted. The result of the model calculation considering the in-media modification for the (a) C and (b) Cu targets, 
together with (c) the fit result with the p — ui interference for the Cu target. The solid lines show the best-fit results. In (a) and 
(b), the shapes of u — > e + e~ (dotted) and p — > e + e~ (dash-dotted) were modified according to the model using the formula 
m(p)/m(0) = 1 — k(p/po) with k = 0.092. In (c), the solid line shows the p — u> interfering shape with the p/u = 0.64 and the 
interference angle = 2.8 radian. The dotted line is the fit result with the typical values of the p/u ratio = 1.0 and the angle = 
2.0 radian Fl. 



of excess was evaluated by subtracting the amplitude for 
the fit function from the data, in the range of 0.6 to 
0.76GcV/c 2 . The obtained ratios of the excess to the u 
peak in the acceptance are estimated to be 0.40±0.04 and 
0.40±0.04 for the C and Cu targets, respectively. The 
obtained yields of each resonance and the excess are given 
in Table [I] This procedure determines the ratio of p to 
u, provided that they follow the relativistic Breit-Wigner 
distribution without any mass modification. After the 
acceptance correction, the 95% C.L. allowed parameter 
regions are obtained as p/uj < 0.15 and p/u < 0.31 for 
C and Cu targets, respectively. The systematic errors 
arises from the uncertainty of the background estimation, 
which amount to 0.09 and 0.21 for the C and Cu targets, 
respectively. The obtained p/u ratios are much smaller 
than unity, as was previously measured in pp interactions 
at the same energy A possible explanation is that 
most of the p are decaying inside the nuclei due to their 
short lifetime; their mass is modified in nuclear media 
and contribute to the excess. 

A comparison has been done for the data with a 
model considering the in-medium mass modification. In 
this model, the mass decreases linearly as a function of 
the density p. in the following relation: m(p)/m(Q) ~ 
1 — k{p/pa) 0, El The parameter k was expected to 
be 16±6% for p and u meson, where po is the nor- 
mal nuclear density 0. The pole mass was modified by 
the above formula according to the density at the decay 
point. The decay-width modification was neglected. A 
Woods-Saxon shape was used for the nuclear density dis- 
tribution; p/ po oc (1 + exp{{r — R)/t))~ 1 , where R = 2.3 
and 4.1 fin, r = 0.57 and 0.50 fm for the C and Cu tar- 
gets, respectively. We assumed that the vector mesons 
were generated at the surface of an incident hemisphere of 



the target nucleus. This assumption is reasonable, since 
we have observed the mass-number dependence of the u 
production cross section as cr(A) oc A 2 / 3 0]. This model 
predicts that the probabilities of p meson decays inside 
a nucleus are 46% and 61% for the C and Cu targets, 
respectively, while those of u are 5% and 9%. 

Based on the model, we modified the shapes of p and u, 
by introducing the shift parameter k which is common for 
the C and Cu targets. We fit again the entire mass region 
of FigFJ using the same procedure as before. The fit 
results, after subtracting the combinatorial background 
and the shapes of ?/ — > e + e _ 7 and u — > e + e~7r° are 
shown in Figl^a) and (b). The spectra for both C and 
Cu targets can be reproduced quite well by this model. 
Confidence ellipsoids for the shift parameter fc, and the 
p/u ratio are also shown in Fig [21 

We obtained the shift parameter of k = 0.092±0.002. 
The best-fit values of the p/u ratio are 0.7±0.1 and 
0.9±0.2 for the C and Cu targets, respectively. Please 
note that here the yield of p and u includes both free 
decays and in-medium decays. It is concluded that the 
observed modification can be understood with the model 
in which the masses of the p and u mesons decrease by 
9% at the normal nuclear density. This value is consistent 
with the theoretical prediction |fij . 

We have neglected in-media mass broadening, which is 
theoretically predicted in and well summa- 

rized in |2.1l2 1| , since it appears that such width broaden- 
ing of p and u, which increases the yield of the high-mass 
tail, does not fit our data [29| . 

It is suggested that the p — u interference is a pos- 
sible explanation of the observed spectral modification 
Such possibility has been inspected by fitting the 
data with the interfering shape, which depends on the 
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FIG. 3: Confidence ellipsoids for the shift parameter k (see 
text), and the production ratio p/ui, for the (a) C and (b) 
Cu targets. The step of each ellipsoid corresponds to Ax 2 =L 
The best fit point is indicated by the cross. 



production ratio p/ui, and the interference angle 9 |25| . 
We scanned the p/ui ratio from 0.2 to 2.6, and 9 from 0.0 
to 3.6 radian. It turned out that the fit was unsuccessful 
for any combination in the scanned regions, and rejected 
at 99.9% C.L. Around the suggested interference angle, 
2 radian the minimum \ 2 was found when 9 was 
2.8 and the p/ui ratio was 0.64. This result is shown in 
FigEfc), together with the curve for 9 = 2.0 radian and 
the p/ui ratio = 1.0 as a reference. 

In summary, we have observed the excess over the 
known hadronic sources on the low-mass side of the ui 
meson peak in the e + e~ invariant mass spectra. This 
demonstrates the spectral-shape modification of p and o> 
mesons in nuclear matter. Without the mass modifica- 
tion of hadrons, the obtained p/uo ratio, which is con- 
sistent with zero, contradicts to the known p/tu ratio at 
this energy. A possible explanation is that most of the 
p mesons decay inside the nucleus, and contribute to the 
excess. The observed excess is understood by the model 
in which the mass of p/ui meson decreases by 9% at the 
normal nuclear density, whereas the p — uj interference 
can not explain the data. 
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